Electron Temperature Evolution in Expanding Ultracold Neutral Plasmas 
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We have used the free expansion of ultracold neutral plasmas as a time-resolved probe of elec- 
tron temperature. A combination of experimental measurements of the ion expansion velocity and 
numerical simulations characterize the crossover from an elastic-collision regime at low initial Fe, 
which is dominated by adiabatic cooling of the electrons, to the regime of high F^ in which inelastic 
processes drastically heat the electrons. We identify the time scales and relative contributions of 
various processes, and experimentally show the importance of radiative decay and disorder-induced 
electron heating for the first time in ultracold neutral plasmas. 
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Ultracold neutral plasmas (UNPs) pj] occupy an ex- 
otic regime of plasma physics in which electron and ion 
temperatures are orders of magnitude colder than in con- 
ventional neutral plasmas. The electron temperature in 
these systems evolves under the influence of many fac- 
tors, which can occur on very different time scales, such 
as disorder-induced heating 0, three-body recombina- 
tion j^], and adiabatic cooling 0, The relative im- 
portance of the various effects depends critically upon 
initial conditions, and this has complicated the experi- 
mental study of the electron temperature Id, 7 , |8|, l o l . llOll 
and lead to much theoretical debate 0, U 5, [Til [ij. 
We present here detailed experimental measurements and 
numerical simulations that untangle the time scales and 
contributions of the various competing effects and char- 
acterize the transition from elastic-collision-dominated to 
inelastic-collision-dominated behavior. 

UNPs are of fundamental interest because they can 
be in or near the strongly coupled regime, which is 
characterized by the existence of spatial correlations be- 
tween particles and a Coulomb coupling parameter F = 

/ (ATTEoakBT) > 1, where T refers to the temperature 
of the particles and a — (47rn/3)^^/'^ is the Wigner-Seitz 
radius. Ions in UNPs equilibrate with F^ ^ 3 [a[il. The 
initial electron temperature is under experimental control 
and can be set such that a naive calculation of Fg sug- 
gests that electrons are also strongly coupled. However, 
electrons rapidly leave the strongly coupled regime due 
to various heating mechanisms [2, M, [13] that are central 
to studies presented here. 

UNPs are created by photoionizing laser-cooled Sr 
atoms 1J| just above the ionization threshold with a 
10 ns laser pulse. The ion temperature is initially a few 
millikelvin, which is similar to the temperature of the 
laser-cooled neutral atoms, but ions heat within one mi- 
crosecond to about 1 K due to disorder-induced heating 
isl 15 1 . The initial electron kinetic energy (Ee) equals 
the difference between the energy of the ionizing photon 
and the ionization threshold. With a tunable pulsed-dye 
laser, 2Ee/ik-Q can be set from 1-1000 K. Electrons ther- 
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malize locally within 100 ns and globally within 1 /is [5[. 
Simple equilibration would set the initial Te = 2£'e/3fcB, 
but we will discuss processes that can change T^- 

The plasma density follows the profile of the neutral 
atom cloud. By adjusting the laser-cooling parameters 
and imaging the cloud in two perpendicular axes, we en- 
sure that the plasma has a spherically symmetric Gaus- 
sian profile, n{r) = noe~^ 1'^'^ . Deviations from spheri- 
cal symmetry, e.g. {px — '^y)/<^x, are less than 5%. Typ- 
ically, the initial no is ~10^^"i~^ 
plasma is quasincutral (n^ ' 

\d — (eofceTe/nee^)^/^ ^ cr, where Ui and rie refer to ion 
and electron density respectively. UNPs are unconfined 
and expand into the surrounding vacuum, and quasincu- 
trality is maintained during the expansion. 

Electron temperature evolution during expansion of 
ultracold plasmas has been studied using various tech- 
niques. Electron plasma oscillations [^ measured the av- 
erage density in order to obtain the rms radial terminal 
velocity of the ions. This showed that essentially all the 
initial electron energy is converted to ion expansion en- 
ergy. In addition, plasmas with lower Eg, and higher no 
(which would imply Fe >1 [3]) resulted in an anoma- 
lously fast expansion. Numerical simulations [3| showed 
that various electron heating mechanisms explained the 
result. However, the relationship between density and 
the electron plasma oscillation used i n [d] has been called 
into question in subsequent work 13, l3|- R^f. [3| probed 
Tf, by measuring the fraction of electrons escaping the 
plasma during a small electric field pulse and inferred 
that for 10 K < 2Ee/3kB < 300 K and no ~ 5xl0"m-3, 
electrons approach a narrow range of electron tempera- 
tures (20 K < Te < 40 K) about 5 /iS after photoioniza- 
tion. Spatially resolved fluorescence detection of the ions 
[^ [^ in a cylindrical plasma measured expansion ener- 
gies similar to what was observed in 6] , but they found 
significant deviations from theoretical predictions that 
perhaps arose because of the lack of spherical symmetry. 
Tonks-Dattner modes, which have resonant frequencies 
that are sensitive to T^, were recently observed in UNPs 
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0. They found agreement between data and a model 
that assumed constant over 40 /is of evolution. So the 
electron temperature evolution in UNPs is obviously a 
complicated problem that remains unsettled. 

The results presented here have several advantages 
over previous studies. We have a spherically symmet- 
ric plasma that allows the application of exact analytic 
results. Doppler broadening of the ion optical absorption 
spectrum [13| provides a calibrated, model-independent 
measure of the ion velocity and overall plasma expan- 
sion with excellent temporal resolution stretching from 
the phase of initial ion acceleration to the onset of termi- 
nal velocity. A combination of experiment and numeri- 
cal simulation allows the contributions of various electron 
heating and cooling mechanisms to be separated as never 
before, and we find excellent agreement between exper- 
iment and theory with no adjustable parameters. We 
also present a systematic study of the full spectrum of 
dynamics observed in current UNP experiments. 

For an experimental probe, the Doppler width of the 
ion absorption spectrum Il3| for the entire plasma mea- 
sures •y/((v~^"Ipy = Vi^rms [iSl, where z is the laser prop- 
agation direction, v is the total ion velocity including 
random thermal motion and expansion (u), the angled 
brackets refer to an average over the plasma density and 
velocity distribution. In a quasineutral Gaussian plasma 
such as a UNP, the electron temperature can be found 
from measurements of Vi^rms due to its sensitivity to u 
and the fact that the expansion acceleration. 



fcs (Te + T,) Vn 
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arises from thermal pressure [3, d, @, 113] ■ 

For high Ef. and low no, which we denote as the "elas- 
tic collisional regime" (Fg < 0.1 ^lij), all collisional pro- 
cesses in a UNP are elastic, the plasma expands adiabat- 
ically, and electrons cool 2l| . This leads to a self-similar 
expansion that preserves the Gaussian phase-space dis- 
tributions and is described by an analytic solution of the 
Vlasov equations [3, Q that was originally derived for 
short-pulse laser experiments [l^j 11^ E^l and applied to 
UNPs in • UNPs provide the first clean realization of 
this analytic solution, and this was shown experimentally 
in [23|. The Vlasov equations do not include a collision 
term, which is appropriate because such a term vanishes 
for a Maxwell-Boltzmann velocity distribution. So in this 
sense, the expansion in the elastic collisional regime can 
also be called "coUisionless" [2o| . 

In this regime, Vi^rms is given by [ISl 



[t^ (Te + r.) + ri^r.]. 
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The characteristic expansion time r^xp is given by Texp = 
v/mjCr(0)2/fcB[T'e(0) + rj(0)], and the electron and ion 
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FIG. 1: Expansion velocity (A) and electron temperature (B) 
for low Fe showing little or no electron heating effects. (A) 
The initial peali densities and sizes are no = 3.5 x 10^^ 
and cr(0) = 1 mm. The self-similar analytic solution (dark 
solid line) provides an excellent fit of the data, with fit Te(0) 
indicated in the legend. The full numerical simulation (dashed 
line) describes the data with no adjustable parameters. (B) 
Electron temperature evolutions (Eq. |3]) are determined from 
parameters of the fits to Vi^rms or by numerical simulation. 
For 2Ee/3kB ~ 25 K, the analytic solution fit and numerical 
simulation for Te{t) deviate significantly, showing the impor- 
tance of inelastic collision effects. 



temperatures follow 



Te/, =T,/,(0)/(l + ^Vdp)• 
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Figure [TJ'\. shows the ion velocity evolution for UNPs 
in the elastic collisional regime. Fits using Eq. [2] take 
Tj/e(0) as fit parameters, while cr(0) is fixed to the value 
found from the images. For higher in Fig. [1] the fit is 
excellent and the extracted values of Te(0) yield 2Ee/3kB 
within experimental uncertainty, confirming that inelas- 
tic collisions are negligible in this regime. The underlying 
electron temperature evolution is shown in Fig. [TJ3. As 
expected, Tg drops due to adiabatic cooling because there 
is no significant electron heating. 

For 2Ee/ikB = 25K in Fig.[Tl the fit TeiO) = 33 K ex- 
ceeds the expected value. This provides evidence that in- 
elastic processes are modifying the electron temperature. 
Many processes are expected to contribute. Within tens 
of nanoseconds, disorder-induced heating (DIH) 0, 25 1 
increases electron kinetic energy by as much as several 
kelvin above E^ in the range of our initial conditions. 
DIH is the conversion of potential energy into kinetic as 
spatial correlations develop, which is an effect of strong 
coupling that has been observed for ions in UNPs 
Three-body recombination (TBR) [s^ populates Rydberg 
levels bound by ~ fceTe and heats the free electrons. 
The total TBR rate varies as Te and can be very 
rapid in ultracold systems [l^. Rydberg-electron colli- 
sions (REG) 1,11 can transfer Rydberg atoms to more 
deeply bound levels and heat electrons further, while ra- 
diative decay (RD) of the Rydberg atoms mitigates this 
heating effect. Direct cooling of electrons through equili- 
bration with ions is negligibly slow for these experiments 
because of the large ion-electron mass difference 27 1. 

To understand the interplay between these various ef- 
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fects, we performed numerical simulations, taking into 
account all relevant heating and cooling mechanisms. 
Our description is based on a particle- in-cell simulation of 
the ions and treats the electrons adiabatically as a fluid 
in a constantly changing equilibrium state 0, DIH 
due to particle correlations is accounted for in the initial 
conditions of both plasma components 0, [l5| , assuming 
a homogeneous To for the electrons and a homogeneous 
Ti for the ions. This treatment neglects the influence 
of correlations on later stages of the plasma dynamics, 
which may be a concern because UNP ions are strongly 
coupled [a , but Q showed that the effects of ion-ion 
correlations on the expansion are negligible. Finally we 
use a Monte-Carlo treatment to describe the formation of 
Rydberg atoms and their subsequent binding energy evo- 
lution, employing known expressions for the rates of TBR 
and REC ^] and RD The TBR rates [3 are well- 

confirmed experimentally at high temperatures, but their 
low temperature validity has been questioned 
as it ultimately has to break down for — > due to its 
strong oc To temperature divergence. By comparing 
experiments and calculations the present study indirectly 
tests TBR theory over a wide range of temperatures. 

The numerical simulation reproduces the expansion 
dynamics in Fig. [1] with no adjustable parameters, and it 
shows that 2iJe/3fcB = 25K is on the border of the "in- 
elastic collisional regime" , which is defined here as initial 
Fe > 0.1 [3]. For this data set, Te{t) stays roughly con- 
stant over the observed evolution time and is not well- 
described by the analytic solution (Eq. [3]). Only TBR 
and REC processes must be included in the simulation 
to accurately describe the expansion. 
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FIG. 2: (A) Expansion velocity and (B) electron tempera- 
ture (B) for moderate initial Fe = 0.2, {2Ee/3kB = 14 K, 
a(0) = 0.9mm, and no(0) = 7 x lO^'^m-^). (A) A fit of the 
expansion to the analytic solution (Eq. (2] solid line) yields 
re(0) = 34 K, which reflects much faster expansion than ex- 
pected for elastic-coUisional dynamics with 2i?e/3fcs = 14K 
(dotted line). Only including TBR and REC in the simula- 
tion (dot-dash line) overestimates the heating. Including RD 
as well (dashed line) brings theory into good agreement with 
experiment. Including DIH increases the initial electron tem- 
perature by about 3 K, but has little effect on the velocity fit. 
(B) The fit Te(0) approximately equals the maximum of the 
actual Te{t) obtained with simulation. 

Fig. displays characteristic plasma dynamics for a 



UNP with lower initial electron energy and higher Fe 
that is in the inelastic collisional regime. The expansion 
is much faster than expected for elastic-collisional-regime 
dynamics with Te(0) = 2Ee/3kB = 14K. Incorporat- 
ing only TBR and REC in the simulation considerably 
overshoots the observed ion velocity evolution because it 
overestimates electron heating. Inclusion of RD, which 
transfers Rydberg atoms to more deeply bound states 
without heating the plasma electrons, produces excellent 
agreement with experiment without adjustable parame- 
ters. 

The analytic solution (Eq. [2]) provides a surprisingly 
accurate description of observed ion velocities, but the 
simulation shows (Fig. [2j3) that the extracted T'e(O) is 
only a phenomenological parameter. Physically, TBR 
and REC heat the electrons in the first ~ 5 /zs, although 
it is mitigated by RD. The increasing Tg, and to a lesser 
extent decreasing density, slows TBR, so that adiabatic 
cooling dominates at later times. The fit Te(0) gives a 
rough estimate of the maximum electron temperature. 
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FIG. 3: (A) Ion velocity and (B) electron temperature for 
Fe ~ 1, {2E^/3kB = 4K, no(0) = 8 x IQI'^ cm-^ o-(O) = 
1.0 mm). (A) The numerical simulation including all heating 
effects matches the data with no adjustable parameters. For 
this relatively high density, the initial heating from DIH slows 
TBR and REC enough to produce an observable effect. The 
analytic solution fit is poor, but yields Te{0) = 49 ± 2K. (B) 
The electron temperature shows drastic heating at early times 
and Te{0) no longer provides a good estimate of Te,max- 



Further decrease of Eg or increase in density pushes 
further into the inelastic collisional regime where a naive 
calculation using = 2Ee/5kB implies Fe >1. In this 
regime, fits using the analytic expansion expression fail 
to reproduce the data (Fig. [3]). Simulations show that 
the approximation of a self-similar Gaussian expansion 
also becomes poor due to the large fraction of ions that 
undergo TBR and the higher rate for this process in the 
higher density central region of the plasma. The mea- 
sured ion expansion velocity indicates the occurrence of 
extreme electron heating from 2Ee/3kB = 4K. For this 
relatively high-density sample DIH makes a significant 
contribution. It quickly raises , which slows recombina- 
tion and leads to a lower T^^max = 53 K. The agreement 
between data and simulation indicates no significant de- 
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viation from classical TBR theory Q . 

Figure [4] summarizes our results for electron heating 
in UNPs. The data is organized according to initial Fg 
[3], and it displays a clear trend in the heating of the 
electrons as previously observed in d, Q . The onset of 
heating occurs at the crossover between the inelastic and 
elastic collisional regime. For initial Fg >1, the expansion 
ceases to be self-similar. 
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FIG. 4: Electron heating summary. The plasma dynamics 
are parameterized well by re(no,-Be). Te.maa; is the maxi- 
mum electron temperature attained during the evolution. The 
heating is negligible in the elastic collisional regime (Fe < 0.1 
[T^). and Te,maa;/(2i7e/3fcs) « 1. In the inelastic collisional 
regime, the heating becomes significant and increases with 
increasing density and decreasing 2Ee/3kB- Beyond Fe ~ 1, 
the expansion ceases to be self-similar. 

We have described the use of spectroscopic measure- 
ments of the ion velocity and numerical simulations as a 
time-resolved probe of the electron dynamics of UNPs. 
It resolves outstanding questions regarding the evolution 
of the electron temperature, and shows that the dynam- 
ics vary greatly for different initial electron kinetic en- 
ergy and plasma density. This work identifies the relative 
contributions and timescales of various processes, demon- 
strates the importance of radiative decay and disorder- 
induced electron heating for the first time in UNPs, and 
shows no discrepancy between observed and theoretical 
rates of TBR. 
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